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>' | 1 This paper is a direct continuation of our publication [3] 

where it was found the exact solution of the Einstein-Maxwell 
equations for two static sources of Reissner-Nordstrom type in the 
state of the physical equilibrium. Here we present the exact solution 
of these equations for the case of two rotating charged sources 
and we proved the existence of the physical equilibrium state also 
for this general case. 

1 Introduction 

In the non-relativistic physics two particles can be in equilibrium if the product 
of their masses is equal to the product of their charges (we use the units G = c = 
1). However, the question on the existence of an analogue of such equilibrium 
state in General Relativity is far to be trivial. Besides the natural mathematical 
complications, in General Relativity arise two different types of the "point" 
centers, namely Black Hole (BH) and Naked Singularity (NS) and one need to 
consider all three configurations BH - BH, NS - NS and BH - NS separately. Yet 
in each case the notion of a physically sensible distance between these objects 
should be defined. 

When the Inverse Scattering Method (ISM) have been adopted for integra- 
tion of the Einstein and Einstein-Maxwell equations it was shown on the exact 
mathematical level that Black Holes and Naked Singularities represent nothing 
else but stationary axially symmetric solitons. Then by the ISM machinery one 
can obtain the infinite families of exact stationary axially symmetric solutions 
of these equations containing such solitons centralized at different points of the 
symmetry axis. The formal construction of such solutions do not represents 
any difficulties apart of the routine calculations in the framework of the well 
developed procedure how to insert a number of solitons into a given background 
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spacetime. However, it is quite intricate task to single out from these families 
the physically reasonable constructions which correspond to a real equilibrium 
states of charged Black Holes and Naked Singularities interacting with each 
other. The point is that in general the stationary axially symmetric solitonic 
solutions possess some features which are unacceptable from the physical point 
of view. These unwanted traits are due the presence in the solutions exotic pe- 
culiarities of the following four types: (i) NUT parameters, (ii) angle deficit at 
the points of the symmetry axis, (iii) closed time-like curves around that parts 
of symmetry axis which are out of the sources and (iv) magnetic charges. The 
NUT parameters are incompatible with asymptotic flatness of the spacetime at 
spatial infinity. The angle deficit is the well known conical singularity violating 
the local Euclidness of space at the points of symmetry axis (it can be treated 
as some singular external strut or string preventing the sources to fall onto or 
to run away each other). Keeping in mind the physical applications we also 
should avoid of any excess of closed timelike curves with respect to those al- 
ready existing inside the sources as an inseparable part of their inner structure. 
Also magnetic charges should be excluded since their presence contradicts the 
Maxwell theory. All four aforementioned phenomena have nothing to do with 
a real equilibrium of the physical bodies and the corresponding equilibrium so- 
lution should be free of such pathologies. To deliberate solution from them one 
need to place the set of the free parameters of the solution under some addi- 
tional restrictions which can be written in the form of some system of algebraic 
equations. The problem is that these equations, even for the simplest case of 
two objects, are extremely complicated and it is difficult to resolve them in an 
exact analytical form in order to see directly whether they have physically ap- 
propriate solutions compatible with the existence of a positive definite distance 
between the sources. 

However, the aforementioned nuisances constitute the real troubles only in 
the general case of the rotating sources. The static case is much more simple 
and it would be not an overstatement to say that for the case of 2 non-rotating 
charged objects the problem have been solved completely. The first indications 
that two static charged masses can stay in real physical equilibrium without 
any struts between them and without any other pathologies came from the 
results of Bonnor [1] and Perry and Cooperstock [2] . In [1] it was analyzed the 
equilibrium condition for a charged test particle in the Reissner- Nordstrom field 
and it was shown that such test body can be at rest in the field of the Reissner- 
Nordstrom source only if they both are either extreme (the charge equal to 
mass) or one of them is of BH type (the charge is less than the mass) and the 
other is of NS type (the charge is grater than the mass). There is no way for 
equilibrium in cases when both masses are either of NS type or both are of BH 
type. The more solid arguments in favour of existence of a static equilibrium 
configuration for the Black Hole - Naked Singularity system was presented in [5] , 
where both sources have been treated exactly, that is no one of the components 
was considered as test particle. These results have been obtained with the 
aid of numerical calculations and three examples of numerical solutions of the 
equilibrium equation have been demonstrated. These solutions can correspond 
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to the equilibrium configurations free of struts, though the authors have not been 
able to show the existence of a positive definite distance between the sources. 
The authors of [5] also reported that a number of numerical experiments for 
two Black Holes and for two Naked Singularities showed the negative outcomes, 
i.e. all tested sets of the parameters was not in power to satisfy the equilibrium 
equation. These findings were in agreement with Bonnor's test particle analysis. 

The explicit analytical resolution of the problem in static case have been 
presented in paper [3] where it was constructed the exact analytic solution of 
the Einstein-Maxwell equations for two sources separated by the well defined 
positive distance and free of struts or of any other unphysical properties 0. We 
showed also that such solution indeed exists only for the BH-NS system and it is 
impossible to have the similar equilibrium state for the pair BH-BH or NS-NS. 
After these results the natural question arises whether the analogous equilibrium 
exists for two rotating sources. It turn out that the answer is affirmative and in 
the present paper we demonstrate the exact equilibrium solution of the Einstein- 
Maxwell equation for two rotating charged objects one of which is a Black Hole 
and another represents a Naked Singularity. 

2 On the general properties of solitonic solu- 
tions 

Metric and electromagnetic potentials for any stationary axisymmetric case in 
cylindrical Weyl coordinates (t, p, z, ip) take the forms: 

2 / 2 2 \ 2 2 / \ 

ds = — / [dp + dz ) + g tt dt + 2g tlp dtd(p + g vv dip , (I) 

9tt9 vv - g 2 tv = -p 2 - (2) 
A t = A t (p,z), A tp = A vt (p,z), A p = 0, A z = 0, (3) 

where all metric coefficients depend only on the variables p, z and the signature 
of the metric corresponds to the following signs of the metric coefficients in 
the diagonal case: / > 0, g vv < 0, g tt > 0. The gravitational solitons (in 
the context of ISM) as exact solutions of pure gravity Einstein equations have 
been introduced in the papers jS] H] • The generalization of this technique for 
the coupled gravitational and electromagnetic fields was constructed in [7]. Its 
more detailed description can be found in [8] and in the book [9]. In this 
generalized approach one starts from some given background solution of the 
Einstein-Maxwell equations and stick into it any desired number of solitons. 
We have to do here with the linear spectral differential equations (Lax pair) for 
the 3x3 matrix function ty(p, z,w), where w is a complex spectral parameter 
independent of coordinates p, z. First of all we choose some background solution 
of the Einstein-Maxwell equations and find from the Lax pair the corresponding 
background spectral matrix ^o(p, z, w). Using the ISM dressing procedure it is 

lr The more details of the derivation of this solution an interested reader can find in paper 
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possible to find explicitly the spectral matrix $F n (p, z, w), corresponding to the 
new solution containing n solitons inserted to the background spacetime and one 
can extract the new metric and new electromagnetic potentials from this ty n . 
The solitonic field added to the background can be characterized by the matrix 
^n^Q 1 — I which is a meromorphic matrix function tending to zero in the limit 
w — > oo and having n simple poles in the complex plane of the parameter w 
(one pole for each soliton). 

In pure gravity case some of these poles can be located at the real axis of the 
ID-plane and the corresponding sources have horizons (that is they are of the BH 
type) while complex poles generate objects with naked singularities. However, in 
the presence of electromagnetic field the formal machinery of the ISM developed 
in [7] in general does not permit for poles to be located at real axis which means 
that by this method one can produce solutions containing sources only of the 
NS typei Nevertheless, also in this case after one obtains the final form of 
solution it is possible to forget the way how it was derived and to continue the 
solution analytically in the space of its parameters in order to get the complete 
family containing solutions with real metric of the physical signature and with 
horizons as well. However, the technical procedure how to do this is simple 
only for the case of one-solitonic solution (that is for the Kerr-Newman case) 
and some simple enough generalization of such procedure was found also for 
two static solitonic objects [TO]. In the general case of two rotating sources (the 
corresponding I2-parametric solitonic solution have been constructed in ([11])) 
this task is much more complicated. Fortunately, there is an effective way to get 
over this difficulty. Because we need to construct solution of the BH-NS type 
we can consider the Kerr-Newman black hole as new background (instead of the 
flat spacetime) and insert to it one soliton of NS type. This is exactly what can 
be done easily with the generating technique proposed in [7] and what we are 
interested in. The exact expressions (in terms of the Ernst potentials) for the 
solution together with the proof that all conditions of the physical equilibrium 
can be satisfied are given below. 

3 The BH-NS system in a stationary state 

Our solution depends on twelve independent real constant parameters 

{m , a , b , g , Mo}, {m s , a s , b s , q s , ^ s }, I = z 2 - Z\ > and c , (4) 

where the parameters with the suffix "0" are related to the background solution 
(black hole) and the parameters with the suffix "s" are the parameters of a 
soliton; all parameters are real. The parameter I (which was chosen positive for 
definiteness) characterizes a coordinate distance between the sources because 
z\ and Z2 determine respectively the location of a black hole and a naked sin- 
gularity on the axis. The constant Co is an arbitrary multiplier in front of the 

2 We said "in general" because it can be shown that the ISM considered in [7J can be adjust 
also to the real uj-poles but only for that special restriction on the parameters of the solutions 
which correspond to the extreme Black Holes. 
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metric coefficient / in ([T} which should be chosen in accordance, e.g., with the 
condition of regularity of the axis at spatial infinity. It is convenient to use two 
functions of the enumerated parameters - a real ao and an pure imaginary a s 
determined by the relations 

ol = ml + bl-al-ql- > 0, a 2 = m 2 s + b 2 s - a\ - q 2 - < 0. (5) 

Our solution is stationary and axisymmetric and depends on two Weyl co- 
ordinates p, z. However, it is more convenient to express it in terms of so called 
"bipolar" coordinates - two pairs of polar coordinate centered respectively at 
the location of a black hole (the coordinates with the suffix 1) and at the location 
of a naked singularity (the coordinates with the suffix 2). Of course, these four 
coordinates should satisfy two additional constraints and each of these four co- 
ordinates can be expressed in terms of Weyl coordinates p, z. The corresponding 
defining relations take the forms 



p = - aoy 1 -y\ = y x l - & 2 y 1 - vl and z = zi+x 1 y 1 =i+zi+x 2 y2, 

(6) 

though, it is worth to note that z\ is not an essential parameter because it 
determines a shift of the whole configuration of the fields and their sources along 
the axis. The inverse relations of bipolar coordinates (xi, y\) corresponding to 
real ctq in terms of Weyl coordinates are 



• Tl = 2 
1 

Za 



' - Zi + cr ) 2 + p 2 + \J{z- z-i- cr ) 2 + p 2 



\J{Z- Z 1 + (Tq) 2 + p 2 - y/{z- Zi - CT ) 2 + p 2 



(7) 



For the coordinates (x2, J/a) corresponding to imaginary a s {a 2 < 0) the similar 
relations are more complicate {z2 =1 + Z\): 




X2 = \l I [{z - z 2 ) 2 +p 2 + a 2 ] + J] [{z - z 2 ) 2 +p 2 + a 2 ) 2 - ail : - : 2 >~ 

in 



■f'2 



(8) 

Sometimes it is convenient also to use instead of pairs of coordinates (x\,y\) 
and (x2,j/2) the pairs of quasi-spherical coordinates (ri,9\) and (r2,02)- 



xi = n - mo , x 2 = r 2 - m s , 

yi = cos 6\ . ?/2 = cos O2 ■ 

The Ernst potentials and metric coefficient / of this solution are: 



(9) 



c — i 2 ( m ° ~ ^0) 2 ( m s ~ ib s ) ^ _ go + i/ip | q s + ip, s ^ 

R\ i?2 Ri R2 
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J_ _ x 2 + ia 8 y 2 + K 1 (x 2 - ayy^) + L 1 (x 1 + a y x ) + S (x 2 + 0^2) 
i?i ~ D 
1 _ si + ia yi + i^ 2 (xi - gpj/i) + L 2 (x 2 + ct s t/ 2 ) 
i? 2 ~ D 



(11) 



D = (xi + ia yi + m - z6 ) [#2 + «a s 2/2 + m 8 - ib s + S (x 2 + cr s y 2 )} - (12) 
- [mo - ib - K 2 (xi - a yi) - L 2 (x 2 + cr s y 2 )} x 
x [m s - ib s - K 1 (x 2 - a s y 2 ) - L 1 (x 1 + <j yi] , 



K ^ = ia s - cr s L ^ = (mp + ib )(m 8 - ib s ) - (q - z/j ) (q s + z/x s ) 

(7o + (7 s + Z' (io - <To)(CTi + <T 2 + Z) 

m - cr r (m - i&o)("i s +ib 8 ) - (qo + i^o) (q b ~ V.) 

-ft 2 = ; 7, ^2 — 



00 + ct s — Z 
So 



(ia s - CT s )((To + cr s — Z) 



cr s (a 2 + aft) (ia s - a s )(a + <7 S - Z) ' 
Y, = (m - ZZ>o) (?s + - ( m « - ib s) {qo + «Mo) , 

Here and in what follows the bar over letters means complex conjugation. 



(14) 
(15) 
(16) 



4 Asymptotics at spatial infinity 

It is easy to see that at spatial infinity, i.e. for p 2 + z 2 — »■ 00, the Ernst potentials 
(10)-(14) have the following asymptotical behaviour: 



s=i- 2{M : iB) + o^), 



* = 9±±*9n. + o(\), 



(17) 

where M and B are the total gravitational mass and total NUT parameter of 
the configuration and Q e and Q m are its total electric and magnetic charges. 
Simple calculations show that these parameters are: 



M = Re 

B = -Im 
Q e = Re 
Q m = Im 



(mp - ib )(l +K 1 +L 1 + go) + (m s - ib s )(l + K 2 + L 2 ) 
\ + S Q -{Kx+Lx){K 2 + L2) 

(m - ib )(l + K 1 +L 1 + Sg) + (m s - ib s )(l + K 2 + L 2 ) 
1 + S - {K! + Li)^ + L 2 ) 

(go + ino) (1 + Ki + £1 + So) + (g s + ^) (1 + K 2 + L 2 ) ~ 
l + So-(ifi+ii)(X 2 + i 2 ) 

(go + ijxp) (l + K 1 +L 1 + Sp) + ( gs + ifj, s ) (l + K 2 + L 2 ) 
l + So- (Ki +L 1 )(K 2 + L 2 ) 



(18) 

(19) 
(20) 
(21) 
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As for the metric coefficient / it should be subject to the constraint / — > 1 in 
the limit r — > oo. This condition can be satisfied by an appropriate choice of 
the coefficient cq. It is easy to see that for this we have to choose 

c = |1 + So - (K x + L X ){K 2 + L 2 )l 2 . (22) 

Then if we would like to have a configuration without NUT parameter and mag- 
netic charge we should impose on the parameters of the solution the restrictions 
B = and Q m = 0. 



5 On the closed time-like curves 

If at some points of the axis (where p = 0) one have gt v ^ 0, this implies (in 
accordance with the relation between the metric coefficients in Weyl coordinates 
9ttgipip — 9ttp — — P 2 ) th & t near these points g vv > 0. Such inequality means that 
near these points of the axis the coordinate lines of the periodic (azimuth angle) 
coordinate, being closed lines, are time-like. To avoid such trouble it is necessary 
to demand that on every part of the axis gt v should vanish. As it follows directly 
from the Einstein - Maxwell equations, on the axis of symmetry p = the value 
fi = gtip/gtt is independent of z and therefore, it is a constant. However, this 
constant can be different on different disconnected parts of the axis separated 
by the sources. Therefore, to exclude the existence of closed time-like curves 
near the axis, first of all we should impose two conditions 

= fli = fl+ , (23) 

where f2_, and 57 + are the constants which are the values of gap/gtt on the 
negative, intermediate and positive parts of the axis respectively. If the condi- 
tions (f2"3"|) are satisfied, the corresponding common value of Q. can be reduced 
to zero by a "global" coordinate transformation of the form t' = t + cup, and 
>p' = (p with an appropriate constant a0 

In order to satisfy the conditions (1231) we need to calculate constants f2 + — 0_ 
and Qi — and put both of them to zero. Calculations show that the first 
constant take simple form: 

0+ - il- = -AB , (24) 



where B is the total NUT parameter given by the formula (|I9[) , while the second 
parameter is much more complicated: 

Cl -n IB "*" x I (1 + 2(?) UMo (25) 

1 (a x +ia s ) + (l-2S)(a x +ia s )W - 



3 It is worth to note here that this transformation actually is not an admissible coordinate 
transformation (because ip is a periodic coordinate but t is not) and its correct interpretation 
is connected with some global restruction of the space-time manifold (some cut-and-past 
procedure) such that the coordinate lines of ip become non-closed while the other closed lines 
appear and they become the coordinate lines of a new asimutal coordinate ip'. 
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The explicit expression for S takes the form: 

s _ <7 (m m s + b b s - q q s - p p s ) 

°oG 2 - ^ - erf - 2a a s ) + (7cr + erg _ ia cr s )(l - a - a s )S 

where So follows from (|14[) and 



(26) 



lo x =m x -ib x +i(a x +ia s ), W a = {I 2 - o% - a 2 ) 2 - 4a 2 a 2 , (27) 
Ho = -2i(Z — cr s — ia + m + i6o)^x + 

(o x + icr s - im x + b x )[(l - a s ) 2 - ctq]+ (28) 
+2(a x + icr s ) [(m + iba)(l - a s + ia ) + <r 2 + a 2 ,] 

X x = (mo + ib )(m x - ib x ) - (q - ifi ) (q* + i^x) ■ (29) 

In these formulas we used the new parameters denoted by the same letters but 
with subscript " x". These new constants are defined by the relations: 

m x = M - m , b x = B - b Q , q x = Q e - q , /i x = Q m - Mo , (30) 
ax + io s = -, — ——= , (31) 

(a s + l(J s ) y/Wo 

T = (mo + ibo) (rn s - ib s ) - (q Q - i^ ) (q s + ifj, a ) - (a s + ia s ) (a - ia s - il) , 
where parameters CQ > M,B,Q e ,Q m have been defined earlier by the relations 



6 On the conical singularity at the axis 

If conditions ([2^]) are satisfied, this does not mean yet that the geometry on 
each part of the axis is regular since at the points of different parts of the axis 
the local Euclidness of spatial geometry still may occur to be violated. This 
behaviour looks like on the surface of a cone near its vortex, where the ratio of 
the length of a circle (surrounding the vortex) to its "radius" is not equal to 
2n. In the solution, on any surface z = const intersecting the axis, the length 
and radius of the circles p = const are represented asymptotically for p — » by 
the expressions L — 2iTy/—g vv and R = \ffp respectively. Using the mentioned 
above relation for metric coefficients in Weyl coordinates gtt9ipip—9ta, = ~P 2 an d 
the condition g tv — > discussed just above, we obtain that the condition of the 
local Euclidness of the geometry at these points, i.e. the condition Lj (2irR) — > 1 
for p — > 0, is equivalent to the following constraints: 

P-=Pi = P+ = l, P = fgtt , (33) 

where P_, Pi and P + are the values of the product fgtt respectively on the 
negative, intermediate and positive parts of the axis. (In accordance with the 
Einstein - Maxwell equations, the product fgtt is constant on a part of the axis 
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where g tlfi — 0, however, these constants again may occur to be different for 
different disconnected parts of the axis.) To obtain conditions which should be 
imposed on the parameters of our solution providing the equations (|33|) to be 
satisfied, we have to analyze the behaviour of metric components on different 
parts of the axis of symmetry. Let's do this. 

Negative semi-infinite part of the axis: {p = 0, — oo < z < z\ — do}. On this 
part of the axis 

x\=z\— z, x 2 = l + zi-z, y 1 =y 2 = —l, (34) 

and the direct calculations of the metric component g u and of the factor / lead 
to the following expressions 

[{z-z.f-alUz-z.-lf-al] 
9tt = — -= , (35) 



f 



c D_D. 



where Co has been defined in (|22]) and D- denotes the value of D (defined by 



(|12p) on the negative semi-infinite part of the axis. From these expressions we 
see that the condition (|33[) is satisfied automatically on the negative semi-infinite 
part of the axis. 

Positive semi-infinite part of the axis: {p = 0, z± + 1 < z < oo}. On this part 
of the axis 

xi=z-zi, x 2 =z-zi-l, yx=y 2 = l, (36) 

and here for the metric component gu and the conformal factor we have the 
similar expressions 

= [(*-*i) 2 -<rg][(s-si-Q 2 -<7 fl 2 ] 
9U c D + D+ ' 1 ^ 

c D+D+ 



1 l( z - Zl Y-al][{ Z - Zl -iy-*l] ' 

where D + denotes the value of D on the positive semi-infinite part of the axis. 
From these expressions we see that the condition (|33p is satisfied identically also 
on the positive semi-infinite part of the axis. 

Intermediate part of the axis: {p = 0, z\ + (Jo < z < z i + t}- On this part of 
the axis 

xi=z—zi, X2=l + zi-z, j/i = l, y 2 = -1, (38) 

and the corresponding expressions for g u and the coefficient / on the axis take 
the forms: 

[{z-z.f-alWiz-z.-lf-^] (l-25\ 2 

gtt - ^m, Itt^J ' (39) 



/= 



c A A 



l( z - Zl y-ol][{z- Zl -iy-ol] 
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where parameter function 5 is the same which have been defined already by 
the formula (|26|) . As follows from these last expression the condition (|33|) on 
the intermediate part of the axis (i.e. the equation P = 1) is equivalent to the 
constraint 5 = 0. 



7 Physical magnetic and electric charges 

To obtain a physical result we have to exclude from the solution magnetic 
charges of both sources. To calculate the physical values of magnetic charges 
(note that formal parameters fio and fi 8 are not physical magnetic charges, they 
coincide with them only in the limit of infinite distance between the sources 
I — > oo) we should consider the magnetic fluxes coming from closed space- like 
surfaces surrounding each charged center and apply the Gauss theorem. In this 
way we can find the physical magnetic charges , ni ph ^ (as well as phys- 
ical electric charges q^ h ^ , qf^ ) of each source calculating the corresponding 
Komar-like integrals. The detailed procedure how to do this have been described 
in the section " Physical parameter of the sources" in paper [3] . The results of 
these calculations are: 

q<* h) = q + Re F , ^ h) = Mo + Im F , (40) 

= qx - Re F , n<? ^ = Mx - Im F , (41) 

where 

p . , (a x + i<T s - im x + b x ) £ Q H (1 + 25) 

*-Wx+*Mx-l 2(a x+l a s ) 2W (a x + ia s ) (1 - 2<5) ' ^ 

and we introduced here the new parameter polynomial: 

Co = [(1 + <J S ) 2 (qx +«^x)+ (43) 

+2 [X x - i(a x + ia s )(l + a s - ia )} (q + i/j, ) . 

The physically acceptable solution corresponds to the restrictions fj,^' 1 ^ = fi{ ph * > = 
0. 



8 Summary for the equilibrium conditions 

Here we present a summary of the conditions which provide an equilibrium 
of two interacting Kerr-Newman sources and which should be satisfied by the 
parameters of the whole configuration. Thus, if we fix an arbitrary constant 
multiplier in the metric coefficient / as it was described above, we have an 
eleven-parameter solution. The first of the equilibrium conditions is the vanish- 
ing of a NUT parameter B. Then from (IH))) we have: 

j I" (m - ibp)(l + K x + Li + go) + (m, - ib s )(l + K 2 + L 2 ) ~\ = 
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The second condition provides the local Euclidness of the geometry on the axis 
(i.e. the absence of the conical points) that is equivalent to the vanishing of the 
parameter <5 which was defined by the expression (|26[) : 

m m s + b Q b s = q Q q s + ^oMs ( 45 ) 

The third condition follows from the absence of the closed time-like curves. 
This demand consists of two restrictions f2+ — 17 _ =0 and f2, — f2_ = 0. Taking 
into account equations fl33} and (gSJ, these two restrictions gives now only one 
independent relation f2j — £1_ =0 which acquire the form (see (|24|) and (12510 : 

W w x S7 x - HoHo (46) 

And finally we have to eliminate the physical magnetic charges H ph \ /if^oi the 
sources. In accordance with (|40)) - (|43l) . bearing in mind that 5 = in agreement 
with (J45j) . we have to add the following two conditions: 



fi + Im 

Mx - Im 



(q x +«Mx) 

(q x + ifj, x ) 



(a x + ia s — im x + b x ) 
2(a x + ia s ) 

(a x + ia s - im x + b x ) 



2Wo(o x +ia s ) 



2(a> 



2W D (a x +icr s ) 



= (47) 
(48) 



9 Evidence of the existing of the physical equi- 
librium 

The easiest way to prove that equilibrium equations have physically acceptable 
solutions is to consider the limit of large distance between the sources. Assuming 
that I is much larger than all other parameters we can try to find an expansion 
for the solution of equations (|44|) - (|48|) with respect to the quantity 1/1. It turns 
out that solution in the form of such expansion indeed exists. To show this we 
assume that six parameters mo, ao, qo, m s , a s , q s are independent of I while the 
four constants bo, fJ-o,b s , Us follows from the equilibrium equations (|33 ]) -(|48 ]) as 
functions of I and of the aforementioned six parameters mo,ao, qo,m s ,a s , q s . 
Then, in the limit of large I, these functions can be expanded with respect to 
I /I. The calculations show that these expansions are: 



m (a s + iy/m^ - a? s - ) 
= ^ 1 L + °\P 



q (a s + iyjm 2 a - a 2 s - q 2 s i , , 

= 1 L + [P 



(49) 

(50) 
(51) 
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q s aa ( 1 



Here, in accordance with the initial inequalities ([5]), m 2 — a 2 — q 2 < then 
all parameters (|4Ul) - (|5li)) are real as it should be. Using this result it is easy 
to obtain from (fT5)) and (|20p the total mass and total electric charge of the 
configuration: 



(V m o ~ a l - 1o - A 



M = m + m s + i ; + O ( - ) , (53) 



h (\/ m - a - it - A ) / 1 



= qo + q s + — ; + O ( ^ ) . (54) 



where 



v/mg - - q 2 (m a q s - m s g ) 



2 



A = v > "° "» v . (55) 

( m o - 9o) ( m s - a ? - q 2 s - V TO ? - a f - 9?) 

We remind that, again in accordance with the inequalities ([S]), the quantity 
TOq — a,Q — <7o i s positive. 

There is one additional important result which is coming from the solution 
of the system (|4"4"l) - (j48|) : the basic (independent of /) parameters m, qo,m s , q s 
must satisfy the relation 

m m s = q q s . (56) 

This equation follows from (|45l) , that is from the absence of the conic syngular- 
ities at the points of symmetry axis. Then the total amount of the equilibrium 
constraints on the initial 11 parameters mo, do, bo, qo,/j,o,m s , a s , b s , q s ,fx 8 ,l 
are (|4"5j) - (|5"2"|) and (|5"6"j) . that is five conditions from five equations (|4"4")) - (|4"5|) as it 
should be. Consequently the solution describing the physical equilibrium state 
of two charged rotating sources contains six arbitrary constants. This is by two 
constants more then the corresponding static solution we found in [3] which is 
natural since in case of rotating sources the two rotation parameters ao and a s 
appeared. 

It is worth to mention that in derivation of the expressions (|49|) -([55 |) we used 
the condition which lead to the essential simplifications of these formulas. 
Another remark we would like to add is the fact that formulas (|4U)) - (j54"f we 
obtained up to the order 1/Z 3 . These additional details showed no surprises and 
there is no of big interest to exhibit them here. 

The last remark relates to the total angular momentum of the configuration. 
To find it it is necessary to expand the Ernst potential £ at spatial infinity 
up to the order 1/r 2 choosing an appropriate spherical coordinates r, 9. Then 
the angular momentum J reveal itself in the summand 2i Jr~ 2 cos 9 in £. We 
calculated it but the expression is rather long and we will not show it here. The 
only interesting fact is that in the limit of large I, discussed above, the total 
angular momentum is 

J = TOo fl o + rnsO-s + O (jj ■ (57) 
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This result together with relations M = mo+m s and Q e = qo+q s following from 
(j53)) - (f54|) in the first non- vanishing approximation and relation tooto s = qo q s 
show that at the infinite separation between the sources they are just usual 
electrically charged rotating Kerr-Newman objects without any NUT annoyance 
and without magnetic charges. 

Of course the analysis presented in this section can be consider as a proof 
of existence of the equilibrium for the distances / from the region l cr < I < oo 
but the critical value l cr can be extracted only from the exact solution of the 
equilibrium equations (l4*4")) - (l4*51) (it is probable that in general l cr is simply zero 
under an appropriate definition of the distance). This task we postpone for a 
future work. 
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